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We study fermionic pairing in an ultracold two-component gas of ^Li atoms by 
observing an energy gap in the radio-frequency excitation spectra. With con- 
trol of the two-body interactions via a Feshbach resonance we demonstrate the 
dependence of the pairing gap on coupling strength, temperature, and Fermi 
energy. The appearance of an energy gap with moderate evaporative cooling 
suggests that our full evaporation brings the strongly interacting system deep 
into a superfluid state. 
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The spectroscopic observation of a pairing gap in the 1950s marked an important experi- 
mental breakthrough in research on superconductivity C/J. The gap measurements provided a 
key to investigate the paired nature of the particles responsible for the frictionless current in met- 
als at very low temperatures. The ground-breaking BCS (Bardeen-Cooper-Schrieffer) theory, 
developed at about the same time, showed that two electrons in the degenerate Fermi sea can be 
coupled by an effectively attractive interaction and form a delocalized, composite particle with 
bosonic character. BCS theory predicted that the gap in the low-temperature limit is propor- 
tional to the critical temperature Tc for the phase transition in agreement with the experimental 
measurements. In general, the physics of superconductivity and superfluidity goes far beyond 
the weak coupling limit of BCS theory. In the limit of strong coupling, paired fermions form lo- 
calized bosons and the system can undergo Bose-Einstein condensation (BEC). The BCS limit 
and the BEC limit are connected by a smooth BCS-BEC crossover, which has been subject of 
great theoretical interest for more than three decades (t2ll3ll?ll5l) . The formation of pairs generally 
represents a key ingredient of superfluidity in fermionic systems and the gap energy is a central 
quantity to characterize pairing regime. 

The rapid progress in experiments with ultracold degenerate Fermi gases @ has opened up a 
unique test ground to study phenomena related to pairing and superfluidity at densities typically 
a billion times below the ones in usual condensed-matter systems. In cold-atom experiments, 
magnetically tuned scattering resonances ("Feshbach resonances") serve as a powerful tool to 
control the two-body coupling strength in the gas (0). Based on such a resonance, a strongly 
interacting degenerate Fermi gas was recently realized (HI). A major breakthrough then followed 
with the creation of Bose-Einstein condensates of molecular dimers composed of fermionic 
atoms (IPl l7Z?ll771li2ll7?l) . which corresponds to the realization of a BEC-type superfluid in the 
strong coupling limit. By variation of the coupling strength, subsequent experiments (|7?ll75l 
176117^ 1771 17^ began to explore the crossover to a BCS-type system. This BEC-BCS crossover 
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is closely linked to the predicted "resonance superfluidity" (I7PIBZ?1[Z7][^ and a "universal" 
behavior of a Fermi gas with resonant interactions (I23ll24|l . The observation of the condensation 
of atom pairs il5l\16l and measurements of collective oscillations il7\\18h support the expected 
superfluidity at presently attainable temperatures in Fermi gases with resonant interactions. 

Our ultracold gas of fermionic ^Li atoms is prepared in a balanced spin-mixture of the 
two lowest sub-states |1) and |2) of the electronic Is^ 2s ground state, employing methods of 
laser cooling and trapping and subsequent evaporative cooling (9). A magnetic field in the 
range between 650 to 950 G is applied for Feshbach tuning via a broad resonance centered at 
Bq ^ 830 G. In this high-field range, the three lowest atomic levels form a triplet of states 
1 2), and |3), essentially differing by the orientation of the nuclear spin (m/ = 1, 0, —1). In the 
resonance region with B < Bq, the s-wave scattering length a for collisions between atoms in 
states |1) and |2) is positive. Here two-body physics supports a weakly bound molecular state 
with a binding energy i?b = /{ma"^), where h is Planck's constant h divided by 2-k and m 
is the atomic mass. Molecules formed in this state can undergo Bose-Einstein condensation 
(0177] 1221 123) • At i? = Bq, the two-body interaction is resonant (a ±oo) corresponding 
to a vanishing binding energy of the molecular state. Beyond the resonance (B > Bq) the 
scattering length is negative (a < 0), which leads to an effective attraction. Here, two-body 
physics does not support a weakly bound molecular level and pairing can only occur due to 
many -body effects. 

Our experimental approach (IPII7?I) facilitates preparation of the quantum gas in various 
regimes with controlled temperature, Fermi energy, and interaction strength. We perform evap- 
orative cooling under conditions (051) where an essentially pure molecular BEC containing 

= 4 X 10^ paired atoms can be created as a starting point for the present experiments. 
The final laser power of the evaporation ramp allows us to vary the temperature T. The Fermi 
energy Ep (Fermi temperature Tp = Ep/ks with Boltzmann's constant /c^) is controlled by 
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a recompression of the gas performed by increasing the trap laser power after the cooling pro- 
cess (l25b . The interaction strength is then varied by slowly changing the magnetic field to the 
desired final value. The adiabatic changes applied to the gas after evaporative cooling proceed 
with conserved entropy (I7?t . Lacking a reliable method to determine the temperature T of a 
deeply degenerate, strongly interacting Fermi gas in a direct way, we characterize the system 
by the temperature T' measured after an isentropic conversion into the BEC limit (E5t . For 
a deeply degenerate Fermi gas, the true temperature T is substantially below our observable 
T' (12(511251) . but a general theory for this relation is not yet available. 

Radio-frequency (RF) spectroscopy was introduced as a powerful tool to study interaction 
effects in ultracold Fermi gases (I27ll28ll29t . Molecular binding energies were measured for ''^K 
atoms (EPb . where also the potential of the method to observe fermionic pairing gap energies 
was pointed out. RF spectroscopy was applied for ^Li atoms to study interaction effects up 
to magnetic fields of 750 G (l2St . One important observation was the absence of mean-field 
shifts in the strongly interacting regime. This effect can be attributed to the fact that, in the 
relevant magnetic-field range, all s-wave scattering processes between ^Li atoms in the states 
1 2) and |3) are simultaneously unitarity limited. This property of ^Li is very favorable for 
RF spectroscopy as it suppresses shifts and broadening by mean-field effects. 

We drive RF transitions from state |2) to the empty state |3) at ~80MHz, and monitor the 
loss of atoms in state 1 2) after weak excitation by a 1 -s RF pulse, using state- selective absorption 
imaging (I7?t . Our experiment is optimized to obtain a resolution of ~ 100 Hz, corresponding 
to an intrinsic sensitivity to interaction effects on the scale of ~5 nK, which is more than two 
orders of magnitude below the typical Fermi temperatures. 

We recorded RF spectra for different degrees of cooling and in various coupling regimes 
(Fig. 1). We realize the molecular regime dX B = 720 G (a = +120 nm). For the resonance 
region, we examined two different magnetic fields because the precise resonance location Bq is 
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not exactly known. Our two values B = 822 G (176b and 837 G il3\\18h may be considered as 
lower and upper bounds for Bq. We also studied the regime beyond the resonance with large 
negative scattering length at B = 875 G (a ^ —600 nm). Spectra taken in a "hot" thermal 
sample atT ^ 6Tp (Tp = 15/iK) show the narrow atomic |2) |3) transition line (upper row 
in Fig. 1) and serve as a frequency reference. We present our spectra as a function of the RF 
offset with respect to the bare atomic transition frequency. 

Spectral signatures of pairing have been theoretically considered i30\\31\U2\\33\\J^ . A clear 
signature of the pairing process is the emergence of a double-peak structure in the spectral 
response as a result of the coexistence of unpaired and paired atoms. The pair-related peak 
is located at a higher frequency than the unpaired-atoms signal as energy is required for pair 
breaking. For understanding the spectra both the homogeneous lineshape of the pair signal 
d31\\33l and the inhomogeneous line broadening due to the density distribution in the harmonic 
trap need to be taken into account S34\i . As an effect of inhomogeneity, fermionic pairing due 
to many-body effects takes place predominantly in the central high-density region of the trap, 
and unpaired atoms mostly populate the outer region of the trap where the density is low (051 
I36ll34t . The spectral component corresponding to the pairs thus shows a large inhomogeneous 
broadening in addition to the homogeneous width of the pair-breaking signal. For the unpaired 
atoms the homogeneous line is narrow and the effects of inhomogeneity and mean-field shifts 
are negligible. These arguments explain why the RF spectra in general show a relatively sharp 
peak for the unpaired atoms together with a broader peak attributed to the pairs. 

We observe clear double-peak structures already at T'/Tp = 0.5, which we obtain with 
moderate evaporative cooling down to a laser power of P = 200 mW (middle row in Fig. 1, 
Tp = 3.4/iK). In the molecular regime (B = 720 G), the sharp atomic peak is well separated 
from the broad dissociation signal (29S, which shows a molecular binding energy of = 
h X 130 kHz = /cb X 6.2 /xK. For B approaching Bq, the peaks begin to overlap. In the resonance 
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region (822 and 837 G), we still observe a relatively narrow atomic peak at the original position 
together with a pair signal. For magnetic fields beyond the resonance, we can resolve the 
double-peak structure for fields up to ~900 G. 

For T'/Tp < 0.2, realized with a deep evaporative cooling ramp down to an optical trap 
power of P = 3.8 mW, we observe a disappearance of the narrow atomic peak in the RF spectra 
(lower row in Fig. 1, Tp = 1.2yuK). This shows that essentially all atoms are paired. In the 
BEC limit (720 G) the dissociation lineshape is identical to the one observed in the trap at 
higher temperature and Fermi energy. Here the localized pairs are molecules with a size much 
smaller than the mean interparticle spacing, and the dissociation signal is independent of the 
density. In the resonance region (822 and 837 G) the pairing signal shows a clear dependence on 
density (Fermi energy), which becomes even more pronounced beyond the resonance (875 G). 
We attribute this to the fact that the size of the pairs becomes comparable to or larger than the 
interparticle spacing. In addition, the narrow width of the pair signal in this regime (lower row 
in Fig. I, B = 875 G) indicates a pair localization in momentum space to well below the Fermi 
momentum hkp = \/2mEp and thus a pair size exceeding the interparticle spacing. 

To quantitatively investigate the crossover from the two-body molecular regime to the fermionic 
many-body regime we measure the pairing energy in a range between 720 G and 905 G. The 
measurements were performed after deep evaporative cooling (T'/Tp < 0.2) for two different 
Fermi temperatures Tp = 1.2 fiK and 3.6 fiK (Fig. 2). As an effective pairing gap we define 
as the frequency difference between the pair-signal maximum and the bare atomic resonance. 
In the BEC limit, the effective pairing gap Az/ simply reflects the molecular binding energy 
i?b (solid line in Fig. 2) (Qjl) . With increasing magnetic field, in the BEC-BCS crossover, Az/ 
shows an increasing deviation from this low-density molecular limit and smoothly evolves into 
a density-dependent many-body regime where hAu < Ep. 

A comparison of the pairing energies at the two different Fermi energies (inset in Fig. 2) 
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provides further insight into the nature of the pairs. In the BEC limit, Az/ is solely determined 
by i?b and thus does not depend on Ep. In the universal regime on resonance, Ep is the only 
energy scale and we indeed observe the effective pairing gap to increase linearly with the 
Fermi energy. We find a corresponding relation hAu ^ 0.2 Ep. Beyond the resonance, where 
the system is expected to change from a resonant to a BCS-type behavior, Au is found to depend 
more strongly on the Fermi energy and the observed gap ratio further increases. We interpret 
this in terms of the increasing BCS character of pairing, for which an exponential dependence 
hAu/Ep oc exp(— 7r/2A;i?|a|) is expected. 

In a further series of measurements (Fig. 3), we applied a controlled heating method to study 
the temperature dependence of the gap in a way which allowed us to keep all other parameters 
constant. After production of a pure molecular BEC (T' < 0.2Tp) in the usual way, we adiabat- 
ically changed the conditions to B = 837 G and Tp = 1.2 /iK. We then increased the trap laser 
power by a factor of nine (Tp increased to 2.5/iK) using exponential ramps of different dura- 
tions. For fast ramps this recompression is non-adiabatic and increases the entropy. By variation 
of the ramp time, we explore a range from our lowest temperatures up to T' /Tp = 0.8. The 
emergence of the gap with decreasing temperature is clearly visible in the RF spectra (Fig. 3). 
The marked increase of Au for decreasing temperature is in good agreement with theoretical 
expectations for the pairing gap energy Q. 

The situation of our experiment was theoretically analyzed for the case of resonant two- 
body interaction The calculated RF spectra are in agreement with our experimental results 
and demonstrate how a double-peak structure emerges as the gas is cooled below T/Tp ^ 0.5 
and how the atomic peak disappears with further decreasing temperature. In particular, the 
work clarifies the role of the "pseudo-gap" regime (|5l l^^ . in which pairs are formed before 
superfluidity is reached. According to the calculated spectra, the atomic peak disappears at 
temperatures well below the critical temperature for the phase-transition to a superfluid. A 



7 



recent theoretical study of the BCS-BEC crossover at finite temperature predicts the phase- 
transition to a superfluid to occur at a temperature that on resonance is only ~30% below the 
point where pair formation sets in. 

We have observed fermionic pairing already after moderate evaporative cooling. With much 
deeper cooling applied, the unpaired atom signal disappeared from our spectra. This observation 
shows that pairing takes place even in the outer region of the trapped gas where the density and 
the local Fermi energy are low. Our results thus strongly suggest that a resonance superfluid 
is formed in the central region of the trap L34) . Together with the observations of resonance 
condensation of fermionic pairs (17 51 li 61) and weak damping of collective excitations il7\\18l . 
our observation of the pairing gap provides a strong case for superfluidity in present experiments 
on resonantly interacting Fermi gas. 
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RF offset (kHz) 

Fig. 1. RF spectra for various magnetic fields and different degrees of evaporative cooling. The 
RF offset (kB X 1 /jK ~ /i x 20.8 kHz) is given relative to the atomic transition |2) |3). The 
molecular limit is realized for i? = 720 G (first column). The resonance regime is studied for 
B — 822 G and 837 G (second and third column). The data at 875 G (fourth column) explore the 
crossover on the BCS side. Upper row, signals of unpaired atoms at T' ^ 6Tf (Tp — 15 /jK); 
middle row, signals for a mixture of unpaired and paired atoms at T' = O.bTp {Tp = 3.4 /xK); 
lower row, signals for paired atoms at T" < 0.2Tp (Tp = 1.2 jiK). Note that the true temperature 
T of the atomic Fermi gas is below the temperature T' which we measure in the BEC limit (see 
text). The solid lines are introduced to guide the eye. 



12 




I ■ 1 ■ 1 1 ■ 

700 750 800 850 900 

Magnetic field (G) 



Fig. 2. Measurements of the effective pairing gap Au as a function of the magnetic field B for 
deep evaporative cooling and two different Fermi temperatures Tp = 1.2/xK (filled symbols) 
and 3.6;uK (open symbols). The solid line shows Au for the low-density limit where it is 
essentially given by the molecular binding energy (E5l) . The inset displays the ratio of the 
effective pairing gaps measured at the two different Fermi energies. 
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RF offset (kHz) 

Fig. 3. RF spectra measured at B = 837 G and Tp = 2.5 /xK for different temperatures T' 
adjusted by controlled heating. The solid lines are fits to guide the eye using a Lorentzian curve 
for the atom peak and a Gaussian curve for the pair signal. The vertical dotted line marks the 
atomic transition and the arrows indicate the effective pairing gap Av. 
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Observation of the Pairing Gap in a Strongly 
Interacting Fermi Gas 
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J. Hecker Denschlag, R. Grimm 

Materials and Methods 

Evaporative cooling, trap frequencies and Fermi energy 

Our optical dipole trap consists of a single laser beam (wavelength 1030 nm) focused to a waist 
of 25//m. Evaporative cooling is performed at a magnetic field of B — 764 G by exponentially 
decreasing the trap laser power P from initially 10.5 W, with a fixed time constant of 460 ms, 
down to a variable final power. About 2x10^ atoms are initially loaded into the trap. After the 
full evaporation ramp down to a final power of 3.8 mW, = 2 x 10^ molecules {N — Axl^f 
atoms) are left in the trap. 

The radial and axial trap frequencies u^, oJz are measured as a function of the laser power P 
and described by 

ujr /2t: = 127 Hz x (P/mW)^/^ 

uJz/2t: = (6015/kG + 0.3P/mW)^/2H2. 

At low laser power, the axial confinement is dominated by the curvature of the magnetic field 
used for Feshbach tuning. 

The Fermi energy Ep — h{3u;^u;zNy/^ for a non-interacting gas is calculated from the trap 
frequencies and the measured number N of paired and unpaired atoms in both internal states. 

Thermometry 

To characterize the entropy of the gas we use the temperature T' that is measured after an adi- 
abatic and reversible (i.e. isentropic) conversion of the gas into the BEC limit. In this regime 
better thermometry is available than in the crossover region. We determine T' at a magnetic 
field of 676 G by fitting the well-known bimodal distribution to in-situ images of the trapped, 
partially condensed molecular cloud [SI]. For the full evaporation ramp, we observe a conden- 
sate fraction of >90%. From this lower limit for the condensate fraction (taking into account 
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interaction effects) we determine an upper limit for the temperature T'/Tbec < 0.4, where the 
critical temperature Tbec for a non-interacting molecular BEC of M = N/2 dimers is given 
by ^bTbec = ^(a;^u;zM/l. 202)^3. Using the relation Tbec = 0.518 for a two-component 
Fermi gas in a harmonic trap, we rewrite the temperature in terms of the Fermi energy and we 
finally obtain T' < 0.2 Tp for our deep evaporative cooling conditions. 

For the relation of the true temperature T in the crossover region to our observable T", a 
general theory is presently not available. However, it was theoretically shown that the isentropic 
conversion of a molecular BEC with T'/Tp <S 1 to a non-interacting Fermi gas leads to a strong 
temperature reduction, following a scaling T/Tp oc (T'/Tp)^ [S2]. In a similar way, we can 
also expect a substantial temperature reduction when the molecular BEC is converted into a 
strongly interacting gas with resonant interactions. 

Experimental details for Figures 1-3 

Fig. 1: The measurements shown in the three rows are taken for different evaporative cool- 
ing and adiabatic recompression conditions. The relevant parameters are summarized in the 
following table: 





final evap. power 


N 


recompression 


Tp 


upper row: no evaporation 
middle row: moderate evap. 
lower row: deep evaporation 


10.5 W (no ramp) 
200 mW 
3.8 mW 


2 X 10*^ 
1.0 X 10^ 
4 X 10^ 


10.5 W (no ramp) 
310mW 
34 mW 


15/iK 
3.4/iK 
1.2/iK 



Fig. 2: All measurements are taken with = 4 x 10^ atoms after full evaporative cooling down 
to a laser power of 3.8 mW (same as in the lower row of Fig 1). The filled symbols refer to the 
case of a subsequent recompression to P = 34mW corresponding to Tp = 1.2fiK, and the open 
symbols refer to a recompression to P = 930mW corresponding to Tp = 3.6yLiK. 



Fig. 3: Full evaporation with same parameters as in the lower row of Fig. 1 and in Fig. 2. After 
evaporation, the trap power is increased to P = 310mW yielding Tp — 2.5//K for A?^ = 4 x 10^ 
atoms. The recompression is performed adiabatically (lower picture) or rapidly with different 
time constants to implement controlled heating. 

RF power: In all measurements, we individually adjust the RF power to obtain a maximum loss 
of r~^40% in state |2). This value is chosen as a compromise between good signal-to-noise ratio 
and minimum perturbation of the system. The RF is generally weak and applied to the sample 
for a long time of 1 s to avoid broadening effects. 
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RF molecular dissociation lineshape in low-density limit 

The radio-frequency z/^p couples the atomic states |2) and |3). For a non-interacting atomic 
gas, the transition frequency z/23 as a function of the magnetic field B is determined by the 
well-known Breit-Rabi formula. In the magnetic-field range of interest the nuclear spin almost 
decouples from the electron spin (Paschen-Back regime) and 1/23 is about 80 MHz. We introduce 
Si'RF = '^RF ~ t^23 as the RF offset from the atomic resonance. 

In the case of a weakly bound molecule, formed with a binding energy Ei, by an atom in 
state |1) and an atom in state |2), the RF can dissociate the molecule if the threshold condition 
/i5z/RF > Eb is fulfilled. Above threshold, the RF-induced dissociation produces two atoms 
with a total kinetic energy of 2Ef;. In the center-of-mass system, where the molecule is at rest, 
energy conservation reads 

hSunF = Eb + 2Ek. 

The lineshape of the continuous RF dissociation signal can be understood in terms of the 
Franck-Condon overlap of the molecular wave function with the wavefunction in the dissoci- 
ation channel (a pair of atoms in the states |1) and |3)). The wavefunction of the molecule 
(dimer of atoms in |1) and |2)) is essentially determined by the scattering length a (or the cor- 
responding binding energy E^ pa /ma^). The dissociation-channel wavefunction in general 
depends on the kinetic energy E^ and the scattering length 013. For a <S \avi\, we find that 
the dependence on 013 has negligible influence on the lineshape. However, for ^Li the dissoci- 
ation channel exhibits a broad Feshbach resonance at ~690 G, which significantly affects our 
dissociation lineshape. A calculation of the lineshape [S3] yields the expression 

f{E) oc E-\E - Ekf/\E -Eh + E')-^ 

where E — h Su^f and E' — /ma\^. 

According to our definition, the effective pairing gap h/S.v in the low-density molecular 
limit is directly given by the maximum of the molecular dissociation signal. From the above 
lineshape f{E) it is straightforward to calculate the signal maximum 

hAu = ^Eb , 

where ^ weakly depends on E' / E^ and varies between the two limits C = 1 and ^ = 4/3 for 
E' -C Eb and E' > Eb, respectively. 
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The magnetic-field dependence of Ai/ that we show as solid line in Fig. 2 is calculated on 
the basis of the above lineshape f{E) and the most recent data for the scattering lengths a and 
ai3 (and the corresponding energies Ei, and E') from the NIST group [S4]. 
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